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^ ■ The ^Si—^Di mixing angle for nucleon-nucleon scattering, ei, is calculated 

(•~^ . to next-to-next-to-leading order in an effective field theory with perturbative 

^^ I pions. Without pions, the low energy theory fits the observed ei well for 

Qv I momenta less than ~ 50 MeV. Including pions perturbatively significantly 

0^ ' improves the agreement with data for momenta up to ~ 150 MeV with one less 

parameter. Furthermore, for these momenta the accuracy of our calculation 

is similar to an effective field theory calculation in which the pion is treated 

O I non-perturbatively. This gives phenomenological support for a perturbative 

^ I treatment of pions in low energy two-nucleon processes. We explain why it is 

necessary to perform spin and isospin traces in d dimensions when regulating 



divergences with dimensional regularization in higher partial wave amplitudes. 
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Effective field tlieory provides a teclinique for describing two-nucleon systems in tlie 
most general way consistent with the symmetries of QCD [|l],0|. In Refs. [0,0, Kaplan, 
Savage, and Wise (KSW) devised a power counting that accounts for the effect of large 
scattering lengths. With this power counting the dimension six four-nucleon operators are 
non-perturbative, while pion exchange and higher dimension operators are perturbative. 
Powers of ap are summed to all orders {p is a typical nucleon momentum, and a is an 
S-wave scattering length). When pions are included in a manner consistent with chiral 
symmetry the expansion is in powers of Q/A where Q = p or 171^^, and A is the range of the 
theory. For p < 171^^/2 (below the pion cut), pions can be integrated out leaving only contact 
interactions. Therefore, the theory without pions is an expansion in powers oi p/mj^. Note 
that for low enough momentum the theory without pions will be more accurate since it is 
not limited by the additional m-j^/A expansion. 

A number of observables have been computed at next-to-leading order (NLO) with the 
KSW power counting. These include nucleon-nucleon phase shifts |0-§[, Coulomb corrections 
to proton-proton scattering , proton-proton fusion [|^ , electromagnetic form factors for the 



deuteron M, deuteron polarizabilities H], np — > ^7 |T0|, Compton deuteron scattering [11 



parity violating deuteron processes [jT2[, and ud —^ vd |T3|. Typically errors are 30%-40% 
at leading order (LO) and of order 10% at NLO indicating Q/A ~ 1/3, or A ~ 400 MeV. 
Since the expansion parameter is fairly large, calculations at next-to-next-to-leading order 
(NNLO) are necessary to achieve accuracy comparable to more conventional approaches. 

In the KSW power counting the leading order diagrams for NN scattering are order 
l/Qi so NNLO corresponds to an order Q calculation. In the theory without pions, several 
of the observables listed above have been computed to NNLO Jl^. In the theory with 



pions the potential pion and local operator contributions to the phase shift in the ^Sq 



channel were calculated at NNLO in Refs. [|15],|T6|. The deuteron quadrupole moment |T^ 
has also been computed at this order. In this paper the '^5'i — ^Di mixing angle, ei, is 
calculated at NNLO in the theory with pions. This calculation provides a clear example of 
an observable for which the theory with perturbative pions does better than the theory with 
only nucleons for momenta of order m^, and without additional parameters. In addition, 
for p ~ mTT the accuracy of this prediction is comparable to a calculation which treats the 
pion nonperturbatively 0. 

The relevant Lagrangian has terms with 0, 1, and 2 nucleons: 



£ = ^Tr (9^S a^St) + ^ Tr KS + m^S^) + N^ [iD^ + ^) N 
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^5-^ ^,t_ ^^Q^t ^tn^AAr ^e5i)^o(35i) , ^2 'Vo^gi) nC'^i), ,T^.,'^?Vo('^i) 



Here qa = 1-25 is the nucleoli axial-vector coupling, S = ^^, / = 131 MeV is the pion 
decay constant, the chiral covariant derivative is D^ = d^ + |(^9^.^^ + ■C^^^O; ^^^ ^ = 
2(^"^q^ + ^^^™g^^)? where rriq = diag{mu,'md) is the quark mass matrix. At the order we are 
working ti;Tr(r7i^) = w{mu + wi^) = m^ = (137 MeV)^. Eq. (|l|) contains two-body nucleon 
operators 

of ^) = (iV^pf ^^)iV)t(iV^/^('^^)iV) + /i.e. , (2) 

where the projection matrices are 

2\/2 4v/^^^T ^ ^ n / J ' ^ ^ 

(i = n -|- 1 is the space-time dimension, and V = v — v . The derivatives in Eqs. (^ and 
(§) should really be chirally covariant, however, only the ordinary derivative is needed for 
the calculation in this paper. Cq , C2 , -D2 > ^^^ C*2 in Eq. (|l]) are normalized so 
that the on-shell Feynman rules in the center of mass frame are 

(jCsi) 

3..X3.. =-.c^^^p^ (') 

q(SD) 

where p is the momentum of the nucleon. From now on the superscript (^5*1) will be dropped. 
Eq. (^ is correct even if spin and isospin traces are performed in n dimensions. 

To regulate ultraviolet divergences it is convenient to use dimensional regularization, 
which respects all the symmetries of the Lagrangian. When using dimensional regularization 
it is necessary to perform spin traces in n dimensions in order not to break rotational 
symmetry. This is important for calculating divergent graphs in higher partial waves. For 
the nucleon theory it is convenient to also continue the isospin traces to n dimensions so 



r( Si) 


=-.cr", 




= -.Bp''mJ, 



that the regulator does not break the Wigner symmetry [jT8| of the lowest order Lagrangian 



|[T9| . Spin and isospin polarization vectors are then normalized so that 



For the scattering NN{ei 
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= d-l = n. 










(5) 


l{ej), i = 


j so calculations 


may 


be 


simplified by 


setting 


e.e* 


n 










(6) 



A more detailed discussion of traces in n dimensions is given in Appendix A. 

To implement the KSW power counting it is useful to use a renormalization scheme 
where the power counting is manifest, such as PDS [|,H or OS pO| , ^ . (In this paper the 



PDS scheme will be used.) In these schemes coefficients of certain four-nucleon operators 
have power law dependence on the renormalization point, /ir, and taking fiR ^ p ^ rrivr ~ 
Q makes the power counting manifest. The size of these coefficients is larger than naive 
dimensional analysis would predict due to the presence of a non-trivial fixed point for a — >■ 
oo. A consequence of this is that bubble graphs with Cq's must be summed to all orders. 
This sums all powers of ap P,|22[]. The '^Si coefficients in Eq. (0) scale as Cq{iir) ~ l/<5, 
C2(yUR)p^ ~ <5°, and m^D2{fiB) ~ Q^. These parameters are fixed by the ^5*1 phase shift at 
NLO. 6*2 is an unknown parameter and enters into the ^Si — ^Di amplitude at order Q. 
This is clear from the beta function for C2 \1^r) in the theory without pions: 

/?r^ = /^.^ c^rv.) = (^)crv.) c^rv.) • (7) 

Solving this equation gives p^ C2 (/^i?) ~ V^ I l^R ~ Q- As discussed below, pions give C^ 
an additional logarithmic dependence on [ir. 
The leading order ^5*1 — ^6*1 amplitude is 

-r 1^ 47r 1 47r , , 

M 7 + «p M60 

This amplitude has a pole at p = i7 corresponding to the deuteron bound state. The 
deuteron has binding energy B = 2.22 MeV, so 7 = yMB = 45.7 MeV. With this boundary 
condition the difference between 7 and the observed scattering length a is obtained from 



perturbative contributions to Co ^0 



Co{^r) = C^'{^r)+cS'\i,r) + ..., (9) 
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FIG. 1. The two order Q^ diagrams that contribute to ei 
the dashed lines are potential pions. 



The solid lines are nucleons and 



where Cq (fJ^RJ ~ Q^- In the PDS scheme the expansion in Eq. @ is necessary to obtain /i^ 
independent amphtudes at each order in Q. This expansion is also necessary to ensure that 



higher order corrections do not give an amplitude with spurious higher order poles p0| . |21 
The S matrix for the ^5*1 and ^Di channels is 2 x 2 and can be parameterized using the 



convention in Ref. p3 



^ = 1 + 



iMp 

2n 
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JSD JDD 
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e^^^'o cos 2ei 



le 



i5o+i52gj^2ei ^ 



e2*^2 cos 2ei 



(10) 



In this parameterization the mixing angle is given by 

Mp A^^ 



sin(2ei) 



2-K 



+ («?)'[^' 



The phase shifts and mixing angle can be expanded in powers of Q/A 



(11) 



51512 



+ 4'^ + ..., 



ei 



+ e-« + ef) + ... 



:i2) 



where the superscript denotes the order in the Q expansion. The phase shifts and mixing 
angles start at one higher order in Q than the amplitudes because of the factor of p in 
Eq. (0). Since A^^ starts at Q^, there is no order Q^ contribution to ei. This is consistent 
with the fact that this angle is much smaller than the ^Si phase shift. In the PDS scheme, 
expressions for ^q ' , §2 \ and e^ were given in Ref. [^. Our main result is the calculation 
of e^ . The NNLO predictions for 5q and 82 are not needed to calculate e^ and will be 



presented in a future publication |2J]. Expanding both sides of Eq. (|TT|) in powers of Q 
givesQ 



^The branch cut for the square root in Eq. ([T3|) is taken to be on the positive real axis. 
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ei is determined by the order Q^ graphs in Fig. |T| and does not involve any free parameters. 
The order Q^ mixing amphtude is [^] 

^5D(o) ^ y2^^(-i)|m.Re[A'(a)] - ^Im[A'(a)]l, 

,,/ X 3 32 ia ./ 1 3 \ , ,^ „. , 

Xia) = + — +i[ — + — ^ ln(l - 2ta) , 



(14) 



where 



a = 



V 



m^ 



(15) 

At order Q, the Feynman diagrams that contribute to the ^Si — ^Di amphtude are shown 
in Fig. ^ In addition to potential pions, at this order the S-wave phase shifts can have 
contributions from diagrams with radiation pions [§]. Performing the energy loop integrals 
using contour integration, potential pions occur when a pole from a nucleon propagator is 
taken. Radiation pion contributions come from taking a pole in a pion propagator. For 
graphs with radiation pions it is necessary to count powers of p ~ Qr = \/M m.„ p5| and 
then scale down to p ~ m.,^. Order Q contributions can come from Q^ and Q^ radiation 



pion graphs [|T^, however these vanish for a ^5*1 — ^Di transition. Soft pion graphs begin at 
order Ql, and for p ~ m.„ are order Q^ p^. Relativistic corrections begin at order Q"^ and 
therefore are not included. 

In dimensional regularization a graph with k loops includes a factor of (/i/j/2)*^*^^~'^^ 
(where the extra 2 is inserted for convenience). Spin and isospin traces will be evaluated in 
d — 1 dimensions for the reasons discussed in Appendix A. Of the graphs in Fig. |^ only e) 
and f ) are divergent in c? = 4 — 2e dimensions. The divergence in f ) is cancelled by a graph 
with the NLO (5"^D2 counterterm^ given by Eq. (5.2) of Ref. |2l|. The p^/e divergence in e) 
is cancelled by the new counterterm 



This is consistent with 5q{p — ;► 0) = tt. The sign of our ^Di state is the opposite of Ref. 
making A^^'^^^ in Eq. (|l^) have the opposite overall sign. 



2 The bare coefficients in Eq. ([l|) are written as 0'°''"' = 5"^C + C*^'^**^. In PDS additional 
finite subtractions are made so that C^"^^^"^ = C^jir) — I]5"C(/iij), see Ref. [^ . 
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FIG. 2. Order Q diagrams for ei. The filled circle is defined in Fig. |^, and the diamonds in b) 
denote insertions of the ^5*1 — ^Si operators with coefficients Cq , C2 or Z?2. 
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TT 



(16) 



10 V87r/2/ V2e 

Note that it is crucial to indicate what constants are subtracted along with the 1/e pole. 
The coupling C2 is determined from a fit to the observed ei. If the extracted value is to 
be used in other calculations, then its exact definition including finite subtractions will be 
needed^. The divergence in Fig. ||e) induces ln(/i/j) dependence in C2 (/iij)- In PDS 



ci^^\f^R) = i^CoifiR) - -j^^o{f^R)yg^p 



3v^^, ,^M^iVj^./4x 



(17) 



where k and A are constants. Note that there is only one unknown in Eq. ( P^ ) since a shift 
in the value of k can be compensated by changing the value of A. 

At order Q the diagrams in Fig. ^ give the following amplitudes in the PDS scheme 



^r\. '-^a ~r *^b ~r *^c ~r *^d ~r *^e ~r *^f 5 



(18) 



where 
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(19) 
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^We have not compared our value of Cj (w^tt) to the value extracted from the deuteron 
quadrupole moment [^ for this reason. 
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The function A* (a) is given in Eq. (Q), and the functions 3^ (a) and 2(a) are given in 
Appendix B. The sum of the amphtudes in Eq. ( piUD is: 

i3 /\/r^2 



1) r 



(20) 
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Mm A^'^h'^ 8 1 

P ln(l-2m)A'(a) + ^3^(a) + A'(a) , 

An J 2 J 

where Ci, C2, and (q are /x/j independent dimensionless combinations of couphng constants: 
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(21) 



Co(/i,j) ' 10 V 87r/2 ; 
Ci and C,2 also appear in the NLO '^Si amphtude (see Eq. (|A^)). (^2 can be ehminated by 



imposing the condition that no spurious double pole should appear in this amplitude 

„2 „2 / s^x 2 



C2^^C.-|4f^Vl0gfl + 



2 - 2/2 14^; ^"^ v^^. 



(22) 



The constant Ci is extracted from a fit to the '^S'l phase shift at NLO. The order Q contribu- 
tion to ei contains one unknown parameter, Cq or C2 {i^r)- This parameter is determined 



by fitting to the value of ei from the Nijmegen partial wave analysis |^6[ at low momentum. 
Results for ei are shown in Fig. |^. The solid line is the Nijmegen result. The order Q 
result in the theory with pions |Q is shown by the dotted line. The result of the order Q"^ 
calculation in the theory with pions is given by the dot-dashed line in Fig. |^. The values 
used in Fig. 11 are 



7 = 45.7MeV 



Ci = 0.2345 , C2 = -0.1038 , Ce = 0.385 



(23) 
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FIG. 3. Predictions for the ^Si — ^Di mixing parameter ei. The sohd Une is the multi-energy 
Nijmegen partial wave analysis |26|. The long and short dashed lines are the order Q^ and Q^ 



predictions in the theory without pions [14|. The dotted line is the order Q prediction in the theory 
with pions from Ref. 0. The dash-dotted line is the order Q^ prediction in the theory with pions. 



The value of (q in Eq. (^3|) corresponds to 



C. 



(5D), 



m. 



-4.56 W 



(24) 



For comparison results have also been shown in Fig. ^ for the theory without pions [[14 



where the prediction for ei begins at order Q^. The long dashed line is the order Q^ result 
and the theory prediction has one free parameter. The short dashed line is the order Q^ 
result which has two free parameters. With one less free parameter, the order Q^ prediction 
of the theory with pions does better than the order Q^ prediction of the theory without pions 
for p > 50 MeV. In fact the theory without pions breaks down around mT^/2, as expected 
since this is where the pion cut begins. It has been noted in the literature |^ that many 
observables may not test the power counting for perturbative pions. As can be seen from 
Fig. ^, the mixing parameter provides an example in which perturbative pions clearly give 
improved agreement with the data. 

The dot-dashed line in Fig. ^ improves over the order Q result for p < 140 MeV. For 
p ~ TTiTr, the error in the order Q^ prediction for ei is ~ 20%. Recall that the mixing angle 



9 



is small and an error of ~ 0.5° is consistent with our expectation for a NNLO calculation. 
It is interesting to ask how sensitive the results in Fig. ^ are to the choice of parameters. 
If we use the ^Si scattering length to fix 7 instead of the deuteron binding energy then the 
order Q^ result (dotted line) increases by ~ 1° for p ~ 171^^. Therefore, the mixing angle is 
quite sensitive to the location of the pole. On the other hand, the NNLO prediction is not 
sensitive to the value of Ci obtained from fitting the '^Si phase shift. This is because e^ in 



Eq. ( |13D depends on the linear combination 

^ = C6-0.56Ci, (25) 

but is insensitive to the orthogonal combination. A change in (i can be compensated by a 
change in (^ while keeping z ~ 0.255. Solutions with the same z give similar predictions, 
for instance, taking ^1 = 0.300 and (q = 0.423 gives an order Q^ phase shift that differs by 
< 0.08° from the one shown in Fig. ^. 

A further test of the convergence of the Q expansion is provided by examining the extent 
to which the amplitude violates unitarity. When Eq. ([TT|) is expanded in powers of Q the 
expression for ei is explicitly real at each order in Q. However, one could insert the NLO 
expression for A^^ and A^^ and the NNLO expressions for A^^ into Eq. (|lT]) and solve 
for ei without making a Q expansion. The resulting ei will have an imaginary part which 
is order Q^ in the power counting. Comparing the imaginary part of ei calculated using 
Eq. dnp to eS^^ + ep^ gives |Im(ei)/(eS^^ + ef^)| < 0.2 for p < 180 MeV, which is of the 
expected size for an order Q^ quantity. Also, for p < rriT^ the ratio l^'^'^'^^^^'^^^'^^l < 0.6, 
which is consistent with an expansion parameter of order 1/2. The agreement of the size of 
these terms with our expectations suggests that the Q expansion is under control. 

In Ref. 0, the mixing angle is calculated using Weinberg's power counting. In this 
approach, momentum power counting is applied to the potential and then the Schroedinger 
equation is solved numerically. Solving the Schroedinger equation with the one pion exchange 
potential is equivalent to summing ladder graphs with potential pion exchange to all orders. 
However, all necessary counterterms are not included, so there is a residual dependence 
on the cutoff. This cutoff dependence can be used to give an estimate of the uncertainty 
in the theoretical prediction due to higher order effects. We will compare our calculation 
with that of Ref. 0, however it is important to keep in mind that Ref. includes graphs 
which are higher order in Q than those in Fig. |^. Ref. also includes A's and more 
parameters are varied in the fit. The results of Ref. are shown in Fig. ^ Varying the 
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FIG. 4. Prediction for ei from Ref. [g]. The fit was done to the partial wave analysis in Ref. |26] 
shown by the solid line. The long dashed line uses the cutoff A = 0.6 nip, the short dashed line 
uses A = TTT-p, and the dotted line uses A = 1.3 mp. 



cutoff between 0.6 rup and 1.3 nip gives an uncertainty of 0.7° at p = ?Ti^. This uncertainty 
is comparable to the error in our fit which differs from the data by 0.5° at p = 171^^. The 
error in our calculation increases for larger values of p because our prediction grows with 
p faster than the observed ei. For these values of p the nonperturbative calculation suffers 
from considerable uncertainty. For a cutoff equal to rrip, the prediction grows with p, but 
with a lower value of the cutoff {0.6 nip) the calculated ei provides better agreement with 
data. It would be interesting to work to one higher order in Q and/or include A's with the 
KSW power counting to see if the agreement with data at higher p improves. At one higher 
order in Q a four derivative four nucleon ^^i — ^Di operator appears. However, using the 
renormalization group its coefficient is determined in terms of Co, C2, and C2 

For momenta p <^ ?7i^, effective range expansions can be constructed for the phase shifts 
and mixing angle. By integrating the pion out of the effective field theory coefficients in 
this expansion can be predicted. In Ref. p8[ coefficients in the expansions of p cotS^ ^°\ 
p cot(fo, and ei are obtained from the order Q^ calculations in Ref. 0]. Ref. p8|] found 
that the effective field theory gives parameter free predictions for the higher coefficients, but 
these did not agree with fits |29| to the partial wave data. However, it is not clear whether 
the extraction of higher order terms in the expansion is accurate enough to test the effective 
field theory PU| . In toy models it has been shown that the convergence of the effective field 
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theory predictions for these coefficients is slow [^. This also seems to be the case when 



the effective field theory is applied to real data. In Ref. ||T6| it was found that the order 
Q corrections to the coefficients of p cot 5'^ "^"^ improve the agreement with the ffi values, 
however the observed convergence is rather slow. 

From the amplitude in Eq. (^) the order Q^ corrections to the momentum expansion of 
ei can be derived. The expansion in the theory without pions takes the form ||T^ 



ei = hi 



po 



^p2 _^ ^2 



+ b2 



pO 



+ ... 



(26) 



where bi and 62 are constants, ei has a cut at p = ±^7, so the momentum expansion of ei 
only converges for p < 7. Clearly it would be more useful to expand a function with better 



analyticity properties. Following Ref. |^ this can be done by parameterizing the S-matrix 

as: 

\ f .... .... \ 

(27) 
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cos ei — sm ei 
sin ei cos ei 



\ / g2i5o 
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^2i52 
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V 



cosei smei 
— sin ei cos ei 



p cot ^0; P^ cot 82, and ei have momentum expansions with radius of convergence mj^/2 rather 
than 7. For low energy expansions these variables should be used. The expressions for (5o,2 
and (fo,2 are the same to order Q. The mixing angle in this parameterization is related to 
the one in Eq. (PUf ) by 

tan(2ei) 2A^'^ 



tan(2e0-^.^^^-^_^^^-^,,_^^^. 
In terms of the amplitudes, the ffist two terms in the Q expansion of ei are 



(28) 
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(29) 



In Fig. ^ we plot the order Q and Q^ effective field theory predictions for ei using the 
parameters in Eq. (|23|). The open circles in Fig. ^ are data from Virginia Tech [^]. The 
stars are the Nijmegen single energy fit to the data PB| whose quoted errors are invisible 



on the scale shown. It seems somewhat strange that the data point at p = 265 MeV from 
Ref. [^ differs from the fit in Ref. [^ by more than eight standard deviations, 
ei has a series expansion in p^: 



ei = gip^ + g2P^ + g?.p^ + ••• • 



(30) 
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FIG. 5. Predictions for the mixing parameter ei defined in Eq. (27). The solid hne is the 



multi-energy Nijmegen partial wave analysis [26|. The dotted line is the NLO prediction in the 
theory with pions from Ref. Q. The dash-dotted line is the NNLO prediction in the theory with 
pions. The open circles are data from Virginia Tech ||3^ and the stars are Nijmegen single energy 



data [26 1 whose quoted errors are invisible on the scale shown. 



Fitting this polynomial to the solid line in Fig. ^ for 7 MeV < p < 50 MeV and weighting low 
momenta more heavily than high momenta gives the values in the first column in Table |. 
To estimate the uncertainty in the extraction of the Qi we varied the range of momentum 
and weighting used in the fit. The value of gi is quite stable, while g2 and g^ varied by 10% 
and 50% respectively. The effective field theory predictions for the coefficients gi are: 
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In each g^ the first term is from the order Q^ diagrams in Fig. |I], while the remaining terms 
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Fit to Nijmegen ei 0{Q^) 0{Q) 



gi{fm^) 0.30 ±0.01 0.55 0.22 

52 (fm^) -2.0 ±0.2 -4.1 -1.5 

53 (fm^) 8.7 ±4.3 28 9.5 

TABLE I. Predictions for the coefficients in a momentum expansion of ei at LO and NLO in 
the effective field theory. 

are from the order Q diagrams in Fig. 0. Using the values in Eq. (^) gives the predictions in 
Table. |. At order Q^ the effective field theory is off by a factor of 2. The order Q corrections 
make the predictions closer to the fit values; the error is ~ 25% for gi and g2, while g^ is 
consistent within error. The effective field theory is converging onto the experimental g^, but 
the errors are somewhat larger than anticipated by the power counting. The convergence 
for terms in the expansion of ei is faster than the convergence in the ^Sq channel. 

To summarize, we have computed the order Q^ correction to the mixing parameter ei. 
The effective theory converges onto the observed ei, and errors are comparable to uncertain- 
ties in alternative approaches where the pion is treated nonperturbatively for p ~ m,r- When 
performing low energy momentum expansions, it is important to use a parameterization of 
the S matrix in which the mixing angle has a convergent expansion for p < mT^/2. The ef- 
fective field theory predictions for the coefficients of this expansion converge towards values 
extracted from a fit to low energy data. In the future, it will be interesting to see if including 
the A or going to one higher order in the Q expansion will provide better agreement for ei 
at p > -m^. 

S.F. was supported in part by NSERC and wishes to thank the Caltech theory group 
for their hospitality. T.M and I.W.S. were supported in part by the Department of Energy 
under grant number DE-FG03-92-ER 40701. 

Note Added in Proof: While this paper was being reviewed, the authors completed a 
NNLO calculation of the phase shifts in the ^S'o, '^S'l, and ^Di channels. Predictions for the 
other P and D wave phase shifts were also examined. We find that in some of these channels 
the KSW expansion exhibits large corrections at NNLO which suggest a breakdown of the 



perturbative treatment of pions. A detailed discussion can be found in the preprint |2^ 
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APPENDIX A: Traces in n Dimensions 



In the standard implementation of dimensional regularization in relativistic theories, the 
spin traces are performed in d dimensions |^4[. For non- relativistic nucleon-nucleon scatter- 



ing the spin traces are often done in 3 dimensions, after which the remaining scalar integrals 
are evaluated in d = n + 1 dimensions. This is in agreement with performing a partial wave 
expansion of the matrix elements using Clebsh-Gordan coefficients; a procedure specific to 
n = 3. This approach provides well-defined results for S-wave transitions. However, when 
higher partial waves are considered it becomes necessary to perform the spin traces in n 
dimensions. To see why consider Fig. ^e), and replace the bubble sum by a single Cq for 
simplicity. The numerator of this graph is proportional to 
pi pi' ^ii' 



p2 ^ 



Tr[aVV™(T* a™ a^ ] g™ g"^ k^ ^ , (Al) 



where k and q are the two loop momenta which run through the pion lines. First consider 
setting (5** /n = 5** /3 in Eq. ( [AID ^"^^ performing the trace in 4 dimensions. At very low 
momentum, the result can be expanded in p/rriT^. When this is done, the amplitude from 
this graph is proportional to a constant for low p. However, for a ^5*1 to ^Di transition 
the amplitude should be proportional to p"^ at low momentum. The constant indicates that 
projection onto ^Si — ^Di was unsuccessful. If we keep the 5** /n in Eq. ([Al|), and perform 
the trace in 3 dimensions then the amplitude is still proportional to a constant for low 
momentum. However, if the trace in Eq. ( |AT| ) is done in n dimensions then the amplitude 
is proportional to p^ as it should be. In the ^5*1 — ^Di calculation the two terms in round 
brackets in Eq. ( |A1| ) have m^/e divergences. These divergences cancel in the difference no 
matter how the expression is evaluated, because there is no operator in this partial wave 
to absorb an m^/e divergence. However, the finite rri^ contributions only cancel when spin 
traces and projection operators are evaluated in n dimensions. Therefore, in this paper all 
spin traces will be performed in n dimensions. 

In Ref. [p!9| , p3| it was pointed out that the nucleon contact interactions with no derivatives 



are invariant under Wigner's SU(4) spin-isospin symmetry for a'- ^°\a^ ^^^ -^ oo. If spin 
traces are performed in n dimensions then it is necessary to treat the isospin traces on the 
same footing, otherwise Wigner symmetry will be broken by the regulator. For this reason, 
isospin traces will also be done in n dimensions. For example, if the order Q^ radiation pion 
calculation in Ref. EH] is performed with spin traces in n dimensions, but isospin traces in 3 
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dimensions then the result is not proportional to l/a(^So) — l/a(^Si). However, in Ref. |T^ 
it was shown that Wigner symmetry implies that the order Qf graphs should be proportional 
to l/a(^So) — l/a{^Si). If all spin and isospin traces are performed in n dimensions then 
the value of individual order Qf graphs changes, but the sum gives the same result as in 
Ref. [g§. 

If the partial wave projection operators are chosen to have the normalization given in 
Eq. d^) then doing the traces in n dimensions does not change any calculations in the theory 
without pions. For S-wave transitions in the theory with pions this convention amounts to 
a change of renormalization scheme, since the difference in evaluating a graph is an overall 
multiplicative factor of the form 1 + 0(e). In PDS, subleading terms in the beta functions 
for coefficients of four nucleon operators are affected. When spin and isospin traces are done 
in n dimensions the NLO ^Si (or ^Sq) amplitude is 



^(0) 



m 



[M-^ 



-11 2 






(7^ - a') 
4a2 



7 



ln(l + 4a2) --tan-^(2a 



a 



(A2) 



where 7 = 'y/rriT^ and d and (2 are given in Eq. (^T]). Different schemes will give different 
expressions for (^1^2, but the amplitude in Eq. ([A2D will remain the same. 



APPENDIX B: Expressions for y and Z 



yia 



In this appendix we give expressions for y and Z which appear in Eqs. (jT^) and (PU]): 

2a 



3 5 



2 3 

(15-4a2), ^^ 
+ ^ „„ , ^ ln(l + a 



ia'^ 



tan ^(a) — 



+i 



80a6 
3 



4a3 ; 

2. _ (3 + 16a^ + 16a^) 
' 32^^7 



+ ^^)tan"^(2a) 



ia"" 



Im 






Li2(- 



4a3 
2a^ + ia 



+ 4a2 



(Bl) 



+ Li2(-2a2 



ta 



to-' 



1 a 
^2^^ 2" 



(3 + 10^2), ,, , , 2A (3 + 16a2 + I6a^) 

^ -^ ln(l + 4a^) + ^ — — ln^(l + 4a^) 



Z{a) 



7 9i 21 



16a;' 

3i 3ia 29a^ 



5 ^ 200 



+ 



40 16a3 40c^2 ■ 49^^ 

3 (16a^ - SOa^ - Aia^ - 15a + 15i) 



H 



128a7 
3a2 9 



15 



16a^ 8a2 



)ln2 



(B2) 



80a5 
-9i + 27a-24ia2 + 78a3-16a^ 



ln(l — ia) 



32a5 
(9 + 48a2 + 48a^)r3 



■ln(l-2za) 



64a6 



, 9 , , , , / 1 + 2ia \ IT 

2 ln2(l - 2za) + 2Li2(-l + 2za) + Li2(-^^^) + - 



2n 



In deriving the formula for Z{a) we found it useful to use reduction formulae due to Tarasov 
35| implemented with the program from Ref. |3^ . 
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